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PREFACE 

One  of  the  most  characteristic  features  of  systems  is  that 
they  undergo  processes.  Hence,  general  systems  theory  oust  concern 
itself  vith  the  theory  of  processes.  It  is  clear  that  a  given 
physical  process  can  he  described  mathematically  in  intrinsically 
different  ways,  each  having  advantages  and  disadvantages. 

The  a.im  of  this  Memorandum  is  to  describe  the  invariant  imbed¬ 
ding  approach  to  the  study  of  processes,  an  approach  which  enables 
us  to  convert  boundary  value  problems  into  initial  value  problems. 
This  is  a  great  computational  advantage  and  is  frequently  of  use 
analytically.  There  are  immediate  applications  in  mathematical 
physics  and  in  the  modem  theory  of  automatic  control,  where  tvo- 
point  boundary  value  problems  for  Euler  equations  abound.  The  Memo¬ 
randum  should  be  of  interest  to  mathematical  physicists,  control 
engineers,  and  operations  analysts. 
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SUMMARY 

Many  of  the  processes  considered  in  modern  physics  and  control 
theory  lead,  from  the  mathematical  viewpoint,  to  nonlinear  two-point 
boundary-value  problems .  These  problems  are  difficult  to  treat 
computationally  and  analytically.  A  goal  of  the  theory  of  invariant 
imbedding  is  to  provide  a  systematic  technique  far  converting  these 
boundary-value  problems  into  initial-value  problems  through  use  of 
appropriate  variables  and  the  employment  of  functional  equation 
techniques . 

A  particle  multiplication  process  is  considered  far  illustra¬ 
tive  purposes.  The  classical  equations  and  the  invariant  imbedding 
equations  are  derived  and  various  interconnections  are  discussed. 

The  paper  is  Intended  to  be  self-contained. 
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I.  PRODUCTION 

Cue  of  the  most  important  features  of  systems  is  that  they 
undergo  processes.  Thus  it  is  that  a  general  theory  of  systems 
must  concern  itself  with  the  theory  of  processes.  The  aim  of  this 
paper  is  to  provide  an  introduction  to  seme  recent  developments  in 
the  theory  of  processes. 

Customarily  when  studying  a  particular  process,  we  focus  our 
attention  on  it  and  attempt  to  discern  those  properties  which  are 
of  interest  to  us.  For  example,  we  may  describe  the  state  of  a 
system  S  by  means  of  a  vector  p,  which  is  viewed  as  a  point  in 
a  phase  space.  The  local  properties  of  the  process  are  given  by 
specifying  that  a  system  in  a  state  p  at  time  n  will  be  in  state 
q  at  time  n  +  1.  This  is  symbolized  by  writing 

(1)  q  =  T(p,n). 

Here  we  are  assuming  that  the  process  is  deterministic  and  that 
cause  and  effect  are  known.  If  the  initial  state  is  c  and  if  the 
state  at  time  n  is  pn,  then  the  process  is  described  by  the  equations 

(2)  P0  =  c, 

(3)  pn+l  =  T(pn'n)'  n-  0,1,2,.... 

We  might,  for  example,  wish  to  determine  whether  the  system  remains 
in  a  certain  region  of  phase  space  over  all  time  or  not. 

One  of  the  key  conceptual  tools  for  studying  such  processes, 
which  are  viewed  as  sequences  of  transformations  of  the  state  of 
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the  system,  is  the  observation  that  the  state  at  time  n  is  a  function 
of  the  initial  state  c  and  the  time  n,  i.e., 

(4)  pn  =  f(c,n). 

Furthermore  this  function  f(c,n)  has  a  semi-group  [1]  properly 
expressed  by  the  equation 

(5)  f(c,m  +  n)  =  f (f (c,n),m). 

This  is  a  symbolic  representation  of  the  fact  that  the  state  at  time 
m  +  n,  the  initial  state  having  been  c,  can  be  viewed  as  having  been 
attained  by  letting  n  units  of  time  pass,  so  that  the  system  is 
transformed  into  the  state  f(c,n),  and  then,  considering  this  as 
a  new  initial  state,  letting  an  additional  m  units  of  time  pass. 
Notice  that  use  of  this  principle  leads  us  to  consider  not  only  the 
single  process  of  duration  m  +  n  and  initial  condition  c,  but  a 
class  of  processes  of  arbitrary  duration  and  arbitrary  initial 
condition. 

In  essence,  invariant  imbedding  is  a  generalization  of  this 
concept.  To  study  a  particular  process  we  imbed  that  process  in 
an  appropriate  class  of  processes.  Then  we  express  the  relationships 
among  the  properties  of  the  members  of  the  class.  These  relations 
then  constitute  an  analytical  description  of  the  processes.  In  this 
way  we  are  frequently  able  to  forge  new  tools  for  their  analytic 
and  computational  study.  As  a  rule  this  leads  to  the  employment 
of  other  variables  such  as  length  or  energy,  in  addition  to  the 
classical  one,  time,  as  the  semi-group  variable. 
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A  discussion  of  the  mathematical  treatment  of  processes, 
including  stochastic,  adaptive,  and  control  aspects  can  "be  found 
in  [2],  A  discussion  of  invariant  Imbedding  with  applications  to 
neutron  transport  theory,  multiple  scattering,  diffusion,  etc., 
is  given  in  [3]  where  references  to  the  work  of  TJeno,  Redheffer, 
Preisendorfer  and  others  vill  be  found.  Basic  references  to  the 
pioneering  efforts  of  S.  Chandrasekhar  can  be  found  in  [4],  Our 
aim  is  to  show  how  boundary -value  problems  can  be  transformed  into 
initial-value  problems,  a  desire  prompted  in  part  by  the  desire  to 
make  it  easier  to  bring  the  power  of  modem  digital  computers  to 
bear  in  significant  classes  of  problems. 

In  subsequent  sections  we  shall  illustrate  these  general 
remarks  by  considering  a  particle  multiplication  process.  The 
physical  background  Is  provided,  the  standard  description  of  the 
process  is  presented,  and  then  the  treatment  via  invariant  imbed¬ 
ding  is  indicated.  Various  interconnections  are  discussed,  and 
some  computational  and  analytical  aspects  are  illuminated.  There 
are  immediate  applications  to  analytical  dynamics  and  automatic 
control  theory.  The  discussion  is  intended  to  be  self-contained. 
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II.  THE  CLASSICAL  TRANSPORT  EQUATIONS  [31 

Let  us  consider  a  particle  process  which  takes  place  within 
a  thin  tube.  The  tube,  assumed  homogeneous ,  extends  along  a  t-axis 
from  t  =  0  to  t  =  T.  At  the  right  end,  t  =  T,  a  steady  stream  of 
c  particles  per  unit  time  is  incident  on  the  tube  and  at  uhe  end 
t  =  0  a  stream  of  w  particles  per  unit  time  is  incident  on  the 
tube.  Within  the  rod  the  particles  interact  with  one  another  and 
with  the  material  of  which  the  rod  is  constructed.  The  flows  within 
the  tube  are  described  by  the  functions  u(t)  and  v(t),  where 

(6)  u(t)  =  the  number  of  particles  per  second  passing 

the  point  t  and  traveling  to  the  right, 

and 

(7)  v(t)  =  the  number  of  particles  per  second  passing 

the  point  t  and  traveling  to  the  left, 

(8)  OstsT, 

The  interactions  which  take  place  are  described  by  assuming  that 
when  a  flow  u  to  the  right  and  a  flow  v  to  the  left  are  incident 
on  a  section  of  tube  of  length  A,  the  increment  in  the  flow  to  the 
right  is  F(u,v)A  +  o(a),  and  the  increment  in  the  flow  to  the  left 
is  G(u,v)A  +  o(A).  Consequently,  we  may  write  the  equations 

(9)  u(t  +  A)  -  u(t)  =  F(u(t),v(t))A  +  o(A), 


(10) 


v(t)  -  v(t  +  A)  =  G(u(t),v(t))A  +  o(A), 
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where,  as  usual,  o(A )  is  a  function  of  A  having  the  property  that 

oo  S5  ^ 

Upon  dividing  Eqs.  (9)  and  (10)  hy  A  and  letting  A  tend  to  zero, 
we  find  the  relations 

(12)  3-  =  F(u,v), 

(1-3)  -  £jr  -  G(u,v). 

Furthermore,  in  view  of  the  assumptions  made  concerning  the  flows 
which  are  incident  on  the  tube,  we  have  the  boundary  conditions 


v(T)  =  c 

(15) 

> 

11 

0 

3 

In  this  manner  we  are  led,  along  classical  lines,  to  study  the  two- 
point  boundary -value  problem  of  Eqs.  (12)  -  (15).  Nonlinear  two- 
point  boundary- value  problems  are  known  to  be  quite  intractable, 
both  analytically  and  c omputat i onally  [51. 
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III.  THE  REFACTION  AND  TRANSMISSION  FUNCTIONS  [6 1 

To  apply  the  idea  of  invariant  imbedding  to  the  particle 
process  discussed  earlier,  ve  consider  the  class  of  processes 
taking  place  in  identical  homogeneous  tubes  of  lengths  T  ;>  0, 
subject  to  incident  fluxes,  at  the  right  end,  of  c  >  0.  The  flux 
incident  on  the  left  end  is  held  fixed  at  a  particular  value,  v. 
The  physical  situation  suggests  that  ve  consider  the  fluxes  which 
emerge  from  the  right  and  left  ends  of  the  various  tubes,  which 
ve  call  the  "reflected"  and  the  "transmitted"  fluxes.  Note  that 
the  reflected  flux  which  emerges  from  the  tube  of  length  T,  due 
to  an  incident  flux  c  on  the  right  end,  is  a  function  of  c  and  T. 
Since  there  is  a  constant  flux  w  which  is  incident  on  the  left 
ends  of  the  tubes,  it  is  true  that  the  fluxes  which  emerge  from 
the  tube  ere  dependent  on  w,  but  we  suppress  this  dependence  since 
it  plays  no  real  role  in  our  considerations.  The  reflected  flux 
is  denoted  by 

(1 6)  r  =  r(c,T). 

Similarly,  ve  denote  the  transmitted  flux  by 

(1?)  t  =  t(c,t). 

First  ve  observe  that  the  fluxes  which  emerge  from  a  rod  of 
zero  length  are  known,  namely, 


(18) 


r(c,0)  =  w, 
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(19)  t(c,0)  =  c. 

Next  we  set  about  relating  the  transmitted  and  reflected  fluxes 
for  rods  of  neighboring  lengths;  mere  precisely,  we  express  the 
flux  reflected  from  rods  of  length  T  +  A  in  terms  of  the  flux  re¬ 
flected  fran  rods  of  length  T.  Oir  basic  relation  is 

(20)  r(c,T  +  A)  =  r(c  +  G(r,c)A,T)  +  F(r,c)A  +  o(A). 

This  equation  expresses  the  fact  that  the  flux  reflected  from  a 
rod  of  length  T  +  A,  due  to  an  incident  flux  c,  consists  of  three 
parts : 

A.  Flux  which  is  reflected  from  a  rod  of  length  T  due  to  an 
input  flux  c  which  is  modified  by  passage  through  a  section  of 
length  A,  and  so  is  of  strength  c  +  G(r,c)A  +  o(A). 

B.  An  increment  in  the  flux  to  the  right  due  to  incident 
fluxes  of  strengths  r  and  c  on  a  tube  of  length  A. 

C.  Fluxes  which  are  proportional  to  the  second  and  higher 
powers  of  A. 

Upon  letting  A  tend  to  zero  in  Eq.  (20),  we  arrive  at  the  desired 
equation 

(21)  rT  =  G(r,c)rc  +  F(r,c), 

where  the  subscripts  denote  partial  derivatives. 

In  general  Eq.  (21)  is  a  first-order  nonlinear  partial  dif¬ 
ferential  equation.  Since  it  is  linear  in  the  partial  derivatives 
r  and  r  ,  it  is  referred  to  as  a  quasilinear  partial  differential 
equation  [71.  It  is  to  be  solved  for  the  reflection  function 
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(22)  r  =  r(c,T), 

subject  to  the  initial  condition  in  Eq.  (3). 

In  a  similar  manner  we  see  that  the  transmission  function 
t(c,T)  satisfies  the  relation 

(23)  t(c,T  +  A)  =  T(c  +  G(r,c)A,T)  +  o(A). 

In  the  limit,  as  A  tends  to  zero,  this  becomes 

(2*0  tt  =  G(r,c)T(;. 

In  addition,  the  transmission  function  T  =  t(c,T)  is  subject  to 
the  initial  condition  in  Eg,.  (19). 

Thus  we  have  succeeded  in  relating  the  reflection  and  trans¬ 
mission  properties  of  tubes  of  neighboring  lengths  and  incident 
fluxes  to  each  other.  In  addition,  the  transmission  and  reflection 
properties  of  a  tube  of  zero  length  are  known.  In  this  manner  we 
may  determine  the  functions  r(c,T)  and  t(c,T)  computationally  for 
a  desired  set  of  values  of  the  variables  c  and  T. 
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IV.  DETERMINATION  OF  THE  IHTERH&L  FLUXES  IN  TERMS 
OF  THE  REFLECTED  AMD  TRANSMITTED  FLUXES 

First  let  us  note  that  if  the  function  r(c,T)  is  known,  then 
the  boundary-value  problem  posed  earlier  is,  in  effect,  transformed 
into  an  initial -value  problem.  For  at  the  right  end  of  the  rod  both 
the  flux  traveling  to  the  left 

(25)  v(T)  =  c 

and  the  flux  traveling  to  the  right 

(26)  u(T)  =  r(c,T) 

are  known.  However,  even  more  may  be  said. 

Consider  an  interior  point  t,  where  the  two  fluxes  are  u(t) 
and  v(t).  We  may  consider  that  a  flux  v(t)  is  incident  on  a  rod 
of  length  t,  giving  rise  to  a  reflected  flux  u(t);  i.e., 

(27)  u(t)  =  r(v(t),t). 

Furthermore,  making  use  of  the  transmission  function,  we  note  that 
under  these  conditions  the  flux  emerging  from  the  left  end  of  the 
tube,  v(0)  =  b,  is 

(28)  b  =  T(v(t),t). 

If  now  the  Eas.  (27)  and  (28)  are  solved  for  u  and  v  in  terms  of 
t,  b,  and  w,  then  we  shall  have  obtained  a  general  solution  of  the 
original  system  of  nonlinear  differential  equations.  Though  these 
observations  are  evident  on  physical  grounds,  we  can  easily  verify 
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that  if  u(t)  and  v(t)  are  tvo  functions  that  satisfy  the  relations 
(27)  and  (28),  where  r(c,T)  and  t(c,T)  are  the  reflection  and  trans¬ 
mission  functions  which  satisfy  Eqs.  (21)  and  (24),  then  the  func¬ 
tions  u(t)  and  v(t)  satisfy  the  original  nonlinear  transport  equa¬ 
tions.  In  more  picturesque  language,  we  may  say  that  a  knowledge 
of  the  reflection  and  transmission  functions  enables  us  to  determine 
the  fluxes  within  the  tube. 

Differentiation  of  Eq.  (28)  with  respect  to  t  yields  the 
relation 


(291  0  ■  V  S  +  Tt 

Upon  comparison  with  Eq.  (24) 


(30) 


dv 

dt 


»  G(u,v), 


we  see  that 


providing  that 


(3D 


Tv  t  0. 


Mext  we  differentiate  Eq.  (27)  with  respect  to  t  which  yields 
/™\  du  dv  , 

<32)  at  '  rv  at  rt 


=  -  r  G  +  r. . 
v  t 

Comparison  with  Eq.  (21)  then  shows  that 
(33)  =  F(u,v). 


Thus  our  assertion  is  established. 
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Let  us  note  that  the  assumption 

(3*0  Ty  d  0 

is  a  reasonable  one,  stating  that  the  transmitted,  flux  is  not  inde¬ 
pendent  of  the  incident  flux.  Furthermore  the  constants  b  and  v 
have  the  physical  interpretations 

(35)  u(0)  =  r(v(o),0)  =  v, 


(36) 


b  =  7(v(0),0)  =  v(0). 
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V.  AN  ILLUSTRATIVE  EXAMPIE 

As  an  illustration  consider  a  neutron  multiplication  process 
in  vhich  a  particle  traversing  a  section  of  a  homogeneous  tube  of 
length  h  has  a  probability  of  ah  +  o(h)  of  being  absorbed  by  the 
medium.  Upon  absorption  the  original  particle  disappears  and  is 
replaced  by  two  daughter  particles,  one  traveling  in  each  direction. 
The  tube  is  of  length  T  and  a  flux  of  average  strength  c  is  incident 
on  the  right  end  of  the  tube  and  no  flux  is  incident  on  the  left 
end.  In  this  case  it  is  easy  to  see  that  the  equations  for  the 
average  internal  fluxes  are  [8] 

(37)  u(t  +  h)  =  u(t)  +  ahv(t)  +  o(h), 

(38)  v(t)  =  v(t  +  h)  +  ahu(t)  +  o(h). 


vhich  becane,  upon  letting  h  tend  to  zero, 
du 


(39) 

(40) 


dt 

dv 

dt 


=  av. 


=  au. 


The  boundary  conditions  are 

(41)  u(0)  =  0,  v(T)  =  c. 

Since 

(42)  F(u,v)  =  av 

and 

(43)  G(u,v)  =  au, 
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the  equations  for  the  reflected  and  transmitted  fluxes  are 

(44)  r^,  =  ac  +  arrc 

and 

(^5)  rT  =  arT(;. 

In  viev  of  the  linearity  of  Eas.  (39)  and  (4o),  and  on  physical 
grounds ,  ve  expect  that  the  reflected  flux  will  he  proportional  to 
the  incident  flux  c, 

(46)  r(c,T)  =  R(T)c. 

If  we  substitute  this  expression  into  Eq.  (44),  we  find 

(47)  §§  c  =  ac  +  aR2c, 

or 

(48)  g  =  a(l  +  R2), 

an  ordinary  differential  equation  for  the  reflection  coefficient 
R(T).  For  an  initial  condition  we  have 

(49)  R(0)  ■=  0. 

The  function  R(T)  can  he  found  explicitly  and  Is 

(50)  R(T)  =  tan  a  T. 

!Hius  for  the  reflection  function  r(c,T)  ve  have 

(51)  r(c,T)  =  c  tan  aT. 


This  is  an  interesting  formula  in  itself  and  shows  that  if  the  length 


of  the  tube  is  sufficiently  great,  jt/(2a)  to  be  exact,  then  the 
reflected  flux  becomes  infinite.  Thus  the  "critical  length"  of  the 
tube  is 


(52) 


T 

crit 


jt 

2a 


In  a  similar  maimer  we  find  that 
(53)  t(c,T)  =  c  sec  aT. 


To  find  the  internal  fluxes  we  use  the  integration  theory 
sketched  above.  First  we  write 


(5*0  u  =  v  tan  at, 

(55)  b  =  v  sec  at. 

Upon  solving  for  u  and  v  we  find 


(56) 

v(t)  =  b  cos  at 

(57) 

u(t)  =  b  sin  at 

the  solutions  of  Eqs.  (39),  (^0),  and  (4l) 
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VI.  RIGOROIE  DERIVATION  OF  THE  EQUATIONS  FOR  THE 
REFLECTED  AND  TRANSMITTED  FLUXES  T91 

In  earlier  sections  we  showed  how  to  deduce  the  equations  for 
the  reflected  and  transmitted  fluxes  using  a  physical  argument.  It 
is  important,  though,  to  have  a  purely  mathematical  procedure.  In 
fact,  as  I.  Bus bridge  points  out  in  her  book  [10],  "The  application 
of  these  (invariance)  principles  is  not  easy  and  until  a  precise 
statement  is  given  of  the  physical  conditions  which  are  sufficient 
to  ensure  their  truth,  any  solution  based  on  these  ought  to  be 
verified  in  another  way."  In  this  section  we  shall  show  how  to 
deduce  the  equations  for  the  reflected  and  transmitted  fluxes  from 
the  equations  for  the  internal  fluxes  in  a  purely  mathematical 
manner.  We  shall  see  that  the  essential  feature  is  the  uniqueness 
of  the  solution  of  the  linear  perturbation  equations  of  the  internal 
flux  equations.  The  basic  idea  behind  the  discussion  is  that  there 
is  essentially  only  one  set  of  perturbation  equations,  regardless 
of  what  is  varied. 

Let  us  consider  the  equations 

(58)  ££  =  F(u,v),  u(0)  =  0, 

(59)  -  =  G(ii,v),  v(T)  =  c. 

In  addition  we  consider  the  same  equations  for  the  functions  U(t) 
and  V(t)  on  the  interval  t  =  0tot  =  T  +  h,  with  the  same  boundary 
condit J  ons , 


(60) 


55  =  f(u,v),  u(o)  =  o, 
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(61)  -  =  G(U,V),  V(T  +  h)  =  c. 

We  nay  refer  these  equations  to  the  interval  t  =  Otot  =  Tby  noting 
that 

(62)  V(T)  =  c  -  V'(T)h  =  c  +  G(U(T),V(T))h 

=  c  +  G(u(T),v(T))h. 

In  the  above  equations  we  have  suppressed  all  terns  involving  powers 
of  h  higher  than  the  first.  For  ease  in  writing  we  shall  continue 
to  do  this  throughout  the  section.  Next  we  introduce  the  perturba¬ 
tion  functions  w(t)  and  x(t)  via  the  relations 

(63)  U(t)  =  u(t)  +  w(t)h, 

(6k)  V(t)  =  v(t)  +  x(t)h, 

which  hold  for  0  *  t  <,  ?. 

If  we  write 

(65)  U(T  +  h)  =  U(T)  +  U'(T)h 

=  u(T)  +  w(T)h  +  F(u(T),v(T))h 


and 

(66)  V(0)  =  v(o)  +  x(0)h, 

we  see  that  we  must  study  v(T)  and  x(0)  in  more  detail. 

We  know  that  the  functions  v(t)  and  x(t)  satisfy  the  linear 


perturbation  equations 
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(67) 

(68) 

As  boundary  conditions  we  have 

(69)  v(0)  =  0 

and 

(70)  x(T)  =  G(u(T),v(T) ). 

The  second  condition  follows  from  Eq.  (62).  In  addition,  we  see 
that  if  we  regard  the  functions  u(t)  and  v(t)  as  functions  of  the 
parameter  c,  then 

<71>  k  UC  =  FU(U'V)UC  +  PV(U’V>V 

(72)  -  k  VC  =  Gu(u’v)uc  +  Gv(u'v)vc* 

Furthermore  the  functions  u£  and  vc  satisfy  the  boundary  conditions 

(73)  uc(0)  =  0, 

(7^)  rc(T)  =  1. 

If  ve  now  assume  that  there  is  a  unique  solution  of  Bp.  (67)  and 
(68),  subject  to  the  conditions  (69)  and  (70)  on  the  interval  [0,T], 
then  we  have  the  desired  results 

(75)  w  =  G(u(T),v(T))uc 

and 


=  Fu(u>v)v  +  ?v(u,v)x, 
-  =  Gu(u'v)w  +  Gv(u>v)x- 
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(76)  x  =  G(u(T),v(T))vc. 

For  T  sufficiently  snail  this  assumption  is  certainly  justified. 
Equation  (65)  now  'becomes 

(77)  U(T  +  h)  =  u(T)  +  ucG(u(T),v(T))h  +  F(u(T),v(T)  )h. 

Upon  writing 

(78)  u(T)  =  r(c,T), 

and  letting  h  tend  to  zero,  we  find 

(79)  35  =  G(r,c)  ^  +  F(r,c). 

This  is  the  desired  equation  for  the  reflection  function.  Equation 
(66)  becomes 

(80)  V(0)  =  v(0)  +  G(u,(T),v(T))vch]t=0  . 

This  leads  to  our  equation  for  the  transmission  function 

(81)  =  G(r,c)  ^  . 

In  the  above  manner  we  see  how  the  mathematical  relations 
between  the  equations  for  the  transmitted  and  reflected  fluxes  and 
those  for  the  internal  fluxes  arise.  These  relationships  may  be 
illuminated  from  a  different  mathematical  viewpoint  as  is  shown  in 
[11],  There,  for  linear  systems,  it  is  shown  that  the  equations 
for  the  internal  fluxes  can  be  considered  as  arising  as  Euler  equa¬ 
tions  for  a  certain  quadratic  variational  problem.  The  equations 
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for  the  reflected  and  transmitted  fluxes  arise  upon  applying  the 
principle  of  optimality  f5]  to  the  same  quadratic  variational 
problem.  Lastly,  let  us  note  that  results  of  the  foregoing  nature 
contain  the  variational  formulas  for  one-dimensional  Green's  func¬ 
tions  and  for  the  characteristic  functions  and  characteristic 
values  ri2"|.  The  same  method  applies  to  time-dependent  and  other 
transport  processes. 


20 


VII.  DISCUSSION 

A  Given  physical  process  can  he  considered  from  various 
mathematical  viewpoints.  It  is  to  be  expected  that  certain  prop¬ 
erties  will  be  more  apparent  from  some  viewpoints  than  from  others. 
Frequently  the  classical  derivations  lead  to  boundary-value  problems 
which  are  difficult  to  treat,  especially  from  the  computational 
viewpoint.  One  way  of  overcoming  these  difficulties  is  through  use 
of  the  functional  equation  techniques  of  invariant  imbedding,  as 
we  have  seen. 

Many  of  the  classical  equations  of  mathematical  physics  T7"l 
as  well  as  those  of  modem  automatic  control  theory  [13]  arise  in 
the  form  of  Euler  equations  associated  with  the  minimization  of 
certain  functionals.  As  a  rule  boundary  conditions  either  are  given 
or  arise  in  the  form  of  free  boundary  conditions.  Thus  it  is  that 
a  major  problem  in  all  of  modem  physics  and  control  theory  is  the 
resolution  of  boundary-value  problems.  Application  of  invariant 
imbedding  to  dynamics  can  be  found  in  [1^,15].  An  alternative  to 
the  usual  Hamilton- Jacobi  theory  is  provided. 

The  considerations  of  this  paper  are  readily  generalized  to 
systems  of  equations  of  order  higher  than  two  [6].  They  are  also 
applicable  to  the  differential-integral  equations  which  arise  in 
neutron  transport  theory  [l6]  and  radiative  transfer  [17]. 

From  the  conceptual  viewpoint  let  us  point  out  that  the 
establishing  of  existence  and  uniqueness  of  solution  for  nonlinear 
two -point  boundary -value  problems  is  difficult.  Invariant  imbedding 
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offers  one  avenue  of  approach  fl8].  The  close  conceptual  relation 
between  invariant  imbedding  and  dynamic  programming  is  evident  f51. 

Finally  let  us  note  that  a  direct  computational  attack  on 
nonlinear  two -point  boundary -value  problems  is  possible  via  quasi - 
linearization  [19]. 
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